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Coherent neutrino-nucleon scattering offers a unique approach in the search for physics beyond 
the Standard Model. When used in conjunction with mono-energetic neutrino sources, the technique 
can be sensitive to the existence of light sterile neutrinos. The ability to utilize such reactions has 
been limited in the past due to the extremely low energy threshold (10 — 50 eV) needed for detection. 
In this paper, we discuss an optimization of cryogenic solid state bolometers that enables reaching 
extremely low kinetic energy thresholds. We investigate the sensitivity of an array of such detectors 
to neutrino oscillations to sterile states. A recent analysis of available reactor data appears to favor 
the existence of such such a sterile neutrino with a mass splitting of |Am s toriio| > 1.5 eV and 
mixing strength of sin 2 2# stcr ii c = 0.17 ± 0.08 at 95% C.L. An array of such low-threshold detectors 
would be able to make a definitive statement as to the validity of the interpretation. 
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INTRODUCTION 

The theory of neutrino oscillations, as described by the 
Pontecorvo-Maki-Nakagawa-Sakata (PMNS) mixing ma- 
trix, provides a simple and well-grounded description of 
the neutrino data obtained thus far [l[. Neutrino oscil- 
lation experiments carried out over the past half-century 
firmly establish the presence of the phenomena to the 
point that the existence of non-zero neutrino masses is 
no longer in question. Current experiments are now en- 
gaged in gaining greater precision of the relevant mixing 
and mass parameters in a manner similar to that which 
was done for the quark sector. 

As the precision of such experiments continues to im- 
prove and the tools by which the data is analyzed increase 
in sophistication, one finds that the emerging picture may 
be more complex than previously realized. For example, 
a recent re-analysis of existing reactor data by Mention 
and collaborators [2[ shows a ~ 3er deviation from the- 
oretical predictions. Though it is possible that the dis- 
crepancy could be due to difficult to calculate Standard 
Model effects, such as weak magnetism [3J, the result also 
appears consistent with the presence of a fourth, sterile 
neutrino. A combined analysis using available reactor 
data, as well as data collected by gallium solar neutrino 
calibration experiments [J, [5[ and the MiniBooNE neu- 
trino data [6| leads to a mass splitting of |A?7i 2 tcrilo | > 1.5 
eV 2 and sin (2<9 S ) 2 = 0.17± 0.08 at 95% C.L. This obser- 
vation appears to be further corroborated by other mea- 
surements, particularly the more recent data collected by 
the MiniBooNE experiment in antineutrino running [7j. 
At this stage it is too early to make any strong claim as to 
the validity of one or all of these observations. Continued 
data collection and scrutiny of systematic uncertainties 
will provide better guidance as to whether new physics 
is at play. 

The most likely beyond the Standard Model explana- 
tion given to the observations made at LSND and Mini- 



BooNE and, more recently, from reanalysis of reactor 
data is the existence of at least one sterile neutrino with 
a mass scale of 1 eV. Ongoing short-baseline measure- 
ments, as well as other complimentary approaches [8j, 
should be able to determine if the data continues to di- 
verge from Standard Model predictions. However, to lay 
claim that the observation is indeed that of a sterile neu- 
trino would almost certainly warrant one or more exper- 
iments with unique signatures to the phenomena. In this 
paper, we propose an alternate approach which makes use 
of oscillometry measurements of neutrino-nucleon coher- 
ent scattering in order to positively confirm or refute the 
existence of sterile neutrinos. 



DETECTION VIA COHERENT SCATTERING 

Coherent scattering offers distinct advantages com- 
pared to other techniques in disentangling the signa- 
ture of sterile neutrinos. First and foremost, coherent 
scattering off nuclei is a neutral current process. Thus, 
any observation of an oscillation structure would indi- 
cate mixing solely to non-active neutrinos. Other meth- 
ods, such as neutrino-electron scattering, must disentan- 
gle the mixing to sterile neutrinos from mixing to active 
neutrinos. The technique becomes even more powerful 
when combined with low energy mono-energetic sources. 
Oscillations to neutrinos at the eV mass scale would man- 
ifest themselves over the length of a few meters (for ~ 1 
MeV neutrino energies). The signature would be quite 
difficult to mimic with typical backgrounds. Finally, the 
cross-section for the process is greatly enhanced thanks 
to the coherent nature of the reaction. 

The use of intense neutrino sources to probe sterile 
neutrinos has been proposed previously in the litera- 
ture [9l4l2l| . The difficulty with all such detection schemes 
is the low energy threshold necessary to detect the sig- 
nature nuclear recoil. Such difficulties are circumvented 



by either resorting to targets with low mass numbcrs- 
considerably lowering the cross-section amplitude and re- 
quiring large mass detectors-or by looking instead at the 
charged current reaction using higher energy neutrinos. 
In this paper, we discuss a low energy threshold detector 
based on cryogenic bolometers that has the capability of 
reaching recoil energy thresholds as low as 10 eV. Such 
detectors re-open the door to neutral current coherent 
scattering as a method for sterile neutrino detection. 

Neutrino-nucleus interactions which are coherent in 
character have the advantage of scaling as A 2 , where A 
is the mass number of the target nucleus. For a target 
nucleus with atomic number Z and neutron number N, 
the cross-section as a function of recoil kinetic energy is 
given by the expression [13f : 



where ao{E v ) = -j^-E^Qyy is the total integrated cross- 
section assuming no energy threshold and f(E v ,To) rep- 
resents the fraction of events above a given threshold en- 
ergy, Tq. In the limit that E v <C Ma, the fraction of 
events above threshold can be written as: 



f(E v ,T ) = (l 
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Any detector hoping to detect such a signal with suffi- 
cient statistics must achieve as low a recoil threshold as 
possible. 

THE 37 Ar SOURCE 
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where Gf is the Fermi coupling constant, Ma is the mass 
of the nucleus, F(q 2 ) is the nuclear form factor, and Q 
is the weak charge, defined by the relation: 



Q w = N-Z(l-4sin 2 e w ) 



Oscillometry-based measurements benefit greatly from 
the use of mono-energetic neutrino sources, since it re- 
duces the measurement to a pure flux-versus-distance 
analysis. Low energy electron capture sources provide the 
most effective and clean source of such neutrinos avail- 
able to date [1J| . A number of such neutrino sources have 
been considered in the literature; a few of them are listed 
in Table HI Historically, two such high intensity source 
(2) have been produced for neutrino studies: a 51 Cr source, 



In our study, we will mainly consider mono-energetic 
electron capture sources, all of which have neutrino en- 
ergies below 1 MeV. The maximum momentum transfer 
for such sources is |g max | < 2E V <C 2 MeV. Since the 
form factor F(q 2 ) —$■ 1 for cases where the scale of the 
momentum probe is much larger than the size of the nu- 
cleus, we can safely ignore this correction factor for our 
analysis. 

The maximum kinetic energy imparted on the nuclear 
recoil depends on the neutrino energy and the mass of 
the recoil target: 
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For a silicon target at 1 MeV, that implies a maximum 
kinetic energy of about 50 eV. For a germanium target 
the maximum kinetic energy would be around 20 eV. 
Such low kinetic energies are why detection of the process 
has been so elusive to date. The fraction of events that is 
detectable by a given experiment depends crucially on the 
inherent threshold of the detector. For a monochromatic 
source of energy E v , the effective cross-section can be 
written as: 
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used by the SAGE and GALLEX experiments 

and an 3 Ar gaseous source used in conjunction with the 

SAGE experiment [j|. 

The 37 Ar source is perhaps the most ideal with respect 
to a future coherent-scattering measurement, for a num- 
ber of reasons: 

• 3 Ar produces a very high-energy, near mono- 
energetic neutrino (90.2% at 811 keV, 9.8% at 813 
keV). 

• With the exception of inner bremsstrahlung pho- 
tons, almost all the energy is carried away by neu- 
trinos, facilitating shielding and enabling the source 
to be extremely compact. 

• Extremely high production yield per reactor target. 

The SAGE collaboration successfully produced such a 
source with a total activity of about 400 kCi to be used 
in conjunction with their gallium solar neutrino detec- 
tor. The source was also very compact, extending 14 cm 
in length and 8 cm in diameter, including shielding 17 1. 
Further reduction in size might be possible, even with in- 
creased activity, making 37 Ar an ideal portable neutrino 
source. 

Despite its clear advantages as a source and its his- 
torical precedent, production of such sources is less 
than ideal. The reaction process by which it is gener- 
ated ( 40 Ca(n,a) 37 Ar) requires a high fast neutron flux 
above 2 MeV, an energy regime where few reactors op- 
erate [18j, [l9( . Production also requires large amounts of 



TABLE I: List of properties of selected electron capture neutrino sources. 



Source 


Half-Life 


Progeny 


Production 


Ey 


Gamma (?) 


37 Ar 
51 Cr 
65 Zn 


35.04 days 

27.70 days 

244 days 


37 C1 

Bly 

65 Cu 


40 Ca(n,a) 37 Ar 
n capture on 50 Cr 
n capture on 64 Zn 


811 keV (90.2%), 813 keV (9.8%) 

747 keV (81.6%), 427 keV (9%), 752 keV (8.5%) 

1343 keV (49.3%), 227 keV (50.7%) 


Inner Brem only 
320 keV 7 
1.1 MeV 7 



CaO and processing in nitric oxide, which makes post- 
production handling difficult. Far less complex to pro- 
duce is 51 Cr, which requires only thermal neutrons cap- 
turing on 50 Cr. However, as a source, the high energy 
gamma produced from the decay of the excited state of 
51 V imposes more shielding requirements. As such, in- 
tense 51 Cr may be less ideal for this investigation, but 
still worth considering given the advantages in produc- 
ing the required activity. 

With it's high energy neutrino emission, 65 Zn is also 
an attractive source for consideration |20|. However, its 
1.1 MeV gamma emission complicates the shielding, so 
this source is not considered further. 



THE DETECTOR 



Historical Precedent 



Bolometric detection of neutrino interactions was one 
of the prime drivers for the development of the low- 
temperature detector community. The idea of search- 
ing for signs of coherent (or "gentle") neutrino scat- 
tering with cryogenic bolometers was first suggested by 
Lubkin 21 1, quickly followed with the first experiment 



design by Niinikoski and Udo [22| for detecting coher- 
ent scattering of neutrinos from an accelerator or reactor 
using 1 cm 3 silicon bolometers at 5 mK with an esti- 
mated energy resolution of 2 eV. Their model did not 
incorporate the heat capacity of the thermometer or the 
thermal coupling of the thermometer to the target, but 
established a low-energy threshold very similar to what 
we propose in this paper. 



The detector requirements for this experiment are ex- 
tremely challenging. Due to the low energy of the neutri- 
nos (< 1 MeV), the recoil energy deposited in the target 
is in the order of tens of eV, while the minimum mass 
needed is hundreds of kilograms. Methods of determin- 
ing the energy deposition from particle interactions in 
a target include measuring the ionization, the scintilla- 
tion, and/or the phonon excitations in the material. For 
nuclear recoils of tens of eV, the fraction of the energy 
deposited by the scattering event that produces free or 
conduction band electrons (the quenching factor) is un- 
known at these energies, and is expected to be very low 
(could be zero for some materials). Thus any readout 
scheme involving ionization channels will be at a severe 
disadvantage. Similar uncertainties hold for the scintil- 
lation yield from nuclear recoils at these energies. An 
additional problem for both ionization and scintillation 
readout is that the energy required to create a single 
electron, electron-hole pair, or scintillation photon from 
a nuclear recoil in most liquid or solid targets is a few eV 
for ionization and tens of eV for scintillation. Thus, even 
if any quanta were produced, Poisson statistics would 
make the measurement of the energy of any given recoil 
event fairly poor. We have therefore focused our atten- 
tion to the measurement of phonons created in the inter- 
action. With mean energies of the order of /xeV, thermal 
phonons provide high statistics at 10 eV and sample the 
full energy of the recoil with no quenching effects. 



Cabrera, Krauss and Wilczek [23[ proposed a multi-ton 
silicon bolometer array using deposited superconducting 
films as a thermometer to detect neutrinos from reac- 
tors and the sun. The desire for large total masses and 
the higher available neutrino energies from these sources 
pushed the detector optimization to an array of kg- 
scalc silicon targets, with thermodynamic-noise-limitcd 
thresholds in the hundreds of eV. 



Rare event searches with cryogenic crystal bolometers 
are being actively pursued by several groups in neutrino 
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physics [24H26J and dark matter searches 
of these experiments have energy thresholds hundreds or 
thousands of times higher than the desired threshold for 
this experiment. This is due to optimizations of the sci- 
ence reach for a given experimental low energy threshold, 
operating temperature, detector mass, and readout tech- 
nology. 

The maximum recoil energy induced in the silicon tar- 
get by our proposed neutrino source is around 50 eV. To 
achieve the required sensitivity, we have taken the ap- 
proach suggested by Niinikoski of using arrays of small 
Si bolometers. We chose, however, to follow the approach 
in [23[ and use transition-edge sensors (TES) as the tem- 
perature readout. We took these ideas and optimized the 
design for this sterile neutrino search. 



Detector Design and Expected Performance 

The threshold for a bolometer is a function of its base- 
line energy resolution. A dimensionless measure of the 
sensitivity of a resistive thermometer at a temperature T 
and resistance R is the quantity a, defined asaE-nJ. 
The energy resolution of a TES bolometer is approxi- 
mately [3lJ 



AS rms = oe 



Uk B T 2 C t ot (3 + 1 



(7) 



where fee is the Boltzmann constant, C to t is the total 
heat capacity of the bolometer, and j3 is the exponent 
of the temperature dependance of the thermal conduc- 
tivity between the bolometer and the refrigerator. To 
unambiguously detect events above the noise from the 
detector, we set the experimental threshold to 7.5 oe- 
For a 10 eV threshold, we then need a detector with 
oe < 1-33 eV, or expressed in terms of the full width at 
half maximum, A^fwhm = 2a/21ii2<7£ < 3.14 eV. 

Assuming a conduction path to the cold bath of the 
refrigerator dominated by Kapitza resistance, /3 = 4, and 
with a temperature T = 15 mK, a 10 eV threshold could 
be attained with a heat capacity C to t < 200 pJ/K. How- 
ever, this model is not complete, as it assumes a per- 
fectly isothermal bolometer. In practice, the various in- 
ternal heat capacity systems of the bolometer are decou- 
pled from each other through internal conductances, and 
thcrmalization times of each separate heat capacity must 
also be taken into account. These internal decouplings 
introduce various sources of noise, degrading the energy 
resolution of the bolometer and consequently requiring a 
smaller heat capacity to attain the desired threshold. 

Fig- ID shows a schematic of the model. The bolometer 
is connected to the cold bath at temperature Tt, through 
a weak thermal conductance G p b- The total heat capac- 



ity can be described by Ct, 



Csi + CrES + C, 



where 



Csi oc T 3 is the theoretical heat capacity of Si given by 
Debye theory, Ctes oc T is the TES heat capacity dom- 
inated by the metal electron system, and C OX cess is the 
heat capacity of impurity bands and two-level systems in 
the crystal. The TES phonon system is assumed to be at 
the same temperature as the silicon phonon system, since 
the sub-micron thickness of the TES makes it incapable 
of sustaining its own thermal phonon population. The 
TES electron system is coupled to the phonon system 
through its electron-phonon coupling conductance G ep . 

Impurities and defects in the Si substrate can lead 
to impurity bands and two-level systems (TLS) in the 
crystal which add to the total heat capacity of the sys- 
tem. The size of this excess heat capacity and its equi- 
libration time with the phonon system depends on the 
specific mechanism involved [32| . Si and Ge crystals 
can be acquired with impurity levels of ~ 10 15 and 
10 10 atoms/cm 3 [33j, respectively. At these levels the 
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FIG. 1: Schematic of bolometer model. The refrigerator acts 
as a cold bath at temperature T},. The Si heat capacity is 
connected by a thermal conductance G p b to the bath. The 
TES is connected by the electron-phonon conductance G cp to 
the Si. A potential excess heat capacity with its coupling are 
shown in dashed outlines. For this study we have assumed 
Goxcgss and/or G eX cess can be made small enough to become 
negligible. 



heat capacity C oxcoss and/or the thermal conductance 
Ge XC ess could be low enough to render them negligible 
for our purposes. Ge clearly holds a large advantage in 
this regard, but its lower maximum recoil energy for co- 
herent neutrino scatters and higher Debye temperature 
make the design of the detector more challenging, thus 
making it our backup if the Debye heat capacity values 
cannot be achieved in Silicon at 15 mK. 

Knaak and Meifiner [34j measured the heat capacity of 
a high-purity 19 g silicon crystal with the following im- 
purity concentrations: B: ~5x 10 12 atoms/cm 3 , O and 
C:~ lx 10 15 atoms/cm 3 . Their measurement was done 
down to 50 mK temperatures, very similar in mass, im- 
purity concentration and temperature to what we wish 
to use in our detector. They measure the time-dependent 
heat capacity on three time scales, initial (< 1 ms), short 
(~ 1 ms), and long (> 0.1 s). The initial time scale 
measures the athermal signal before thermalization and 
is less relevant to our present discussion. For the short 
time constant, they report a measured heat capacity con- 
sistent with the Debye prediction. For the long time con- 
stant, they see a larger heat capacity which they attribute 
to the heat capacity of the gold, carbon, and epoxy used 
to instrument their sample. It is possible that some of 
the excess heat capacity in their measurement was not 
due to the instrumentation attached to the sample. The 
time constant of our proposed detector design is 50 ms, 
in between these two reported measurement scales. More 
measurements should be done at these low temperatures 
to study the heat capacity dependence of Si with different 
species and concentrations of impurities. 

The Knaak and Meifiner results suggest that low heat 



capacities are achievable. In our design, we benefit from 
the fact that our TES thermometer makes up about half 
of the total heat capacity, which allows us to tolerate 
some excess heat capacity from the silicon target. We 
also have the option to optimize using less mass per 
bolometer, and trading off overall experiment mass for 
lower threshold per detector. 

For this study, we will assume the excess heat capac- 
ity is negligible, and optimize the heat capacity of the 
thermometer Gtes- the electron-phonon coupling G ep , 
and the Si heat capacity Gs; to obtain the desired 10 eV 
threshold with the highest possible target mass. We make 
the following assumptions: 

• Each detector is a Si cube ranging in mass from 
20-100 g. The heat capacity is determined from 
Debye theory. 

• The conductance between the Si and the cold bath, 
G p b; can be engineered to give a desired value. The 
value is chosen to give a thermal impulse response 
time of 50 ms as measured by the thermometer 
readout. 

• The thermometer is a Mo/Au TES bilayer with 
a superconducting transition engineered to a spe- 
cific temperature between 10-100 mK. Mo/Au 
TES X-ray detectors have achieved resolutions of 
A£ F whm = 2 eV [33. 

• The TES heat capacity and electron-phonon cou- 
pling are taken from the literature and are a func- 
tion of the chosen volume of the TES and the tem- 
perature. 

Given these general assumptions, several combinations 
of detector mass and transition temperature were tested 
for both Si and Ge targets, scaling the TES volume to 
obtain the best e nerg y resolution, following the theoreti- 
cal framework of [36[ . The TES volume is a compromise 
between two competing interests: having a small TES 
heat capacity, and having a fast thermal link between 
the TES and the Si or Ge target. The optimum vol- 
ume corresponds to a TES heat capacity that is roughly 
equal to the target. An important quantity is the ratio 
Gep/Gpb- As long as this ratio is > 100, the TES remains 
in quasi-thermal equilibrium with the target throughout 
a pulse (except for the initial athcrmal phase on the order 
of 1 ms). If needed, one can make G p b smaller and gain 
energy resolution (and lower threshold) at the expense of 
slower signals. 

The results of our models are shown in Fig. [2j The 
plotted threshold is calculated as 7.5 oe- Due to practi- 
cal limitations in refrigeration and considering the read- 
out necessary for the size of the experiment, we focus on 
a transition temperature of 15 mK, with the refrigera- 
tor base temperature at 7.5 mK. At this temperature, 
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FIG. 2: Calculation of threshold for Si and Ge targets from 
20-100 g at different operating temperatures. The lowest line 
for each target material is the 20 g line. The model (see Fig[TJ 
takes into account the heat capacity of the TES and the tar- 
get, the internal thermal fluctuation noise between the target 
and the TES thermometer, the electronics noise, the John- 
son noise from the TES and its bias resistor, and the phonon 
noise between the target to the bath. The volume of the TES 
was scaled to give the best energy resolution at 15 mK. The 
horizontal dashed line marks the desired 10 eV threshold, cor- 
responding to an energy resolution Aj^fwhm = 3.14 eV. The 
vertical dashed line marks the desired operating temperature 
of 15 mK. For Si, a 50 g target meets the requirements. For 
Ge, a 20 g target meets the requirements. 



a low-energy threshold of 10 eV can be obtained with 
bolometers with 50 g of Si or 20 g of Ge. A 50 g Si 
target sees about twice the rate of neutrino coherent re- 
coil events as a 20 g Ge target when both have a 10 eV 
threshold. We will thus focus on Si. Model parameters 
for the Si detector are given in Table [TTJ 

The natural decay time of the bolometer Gtot/G p b = 
436 ms. Electro-thermal feedback J3l| from the TES 
speeds up the response time to roughly 50 ms. The TES 
readout is bandwidth-limited by an inductor which criti- 
cally damps the system, causing a further speedup in the 
response. The decay time of recoil events is reduced from 
the 50 ms decay with no inductor to a 30 ms decay with 
the inductor. Fig. @]shows a simulation of 10-50 eV neu- 
trino coherent scatters in a 50 g Si bolometer. The pulses 
arc clearly separated from the noise, and the energies of 
the different events are clearly separated by eye. 

The total heat capacity is on the order of 460 keV/mK, 
and given a transition width of around 1 mK for a TES, 
we estimate that the bolometer will have fairly linear re- 
sponse up to hundreds of keV. Higher energies will have 
a non-linear response but will retain significant energy 
resolution. This large energy bandwidth will help un- 
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FIG. 3: Zoom-out of Fig. [2] Our desired operating point is 
designated by the vertical and horizontal dashed lines. The 
legend for both plots are the same. The ovals designate 
the region of operation for the following low-temperature ex- 
periments: CUORE [H, EDELWEISS [H, CRESST [37|, 
CDMS [28fl , X-ray microcalorimeters for astrophysics [35|] , and 
MARE (proposed) [2g|. Note that these experiments use dif- 
ferent combinations of thermal or athermal measurements, 
TES or NTD thermometry, and mass per detector unit, so 
cannot be directly compared to the design presented in this 
paper. 



derstand the background in our experiments, and enable 
other rare event searches such as limits on the neutrino 
magnetic moment and dark matter interactions. 
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FIG. 4: Simulated current readout for 10-50 eV recoils us- 
ing the model parameters in Table [IT] The current has been 
multiplied by -1 to make the pulses positive. Noise sources 
modeled are: the phonon noise between the target to the bath, 
the internal thermal fluctuation noise between the target and 
the TES thermometer, the Johnson noise from the TES and 
its bias resistor, and the electronics noise. The modeled 10 eV 
pulses are clearly separated from the noise. 



TABLE II: Model parameters for a 50 g Si target coupled to 
a Mo/Au TES operated at 15 mK. The Si target is a 28 mm 
cube, and the TES is an 25 mm x 2 mm film 600 nm thick 
deposited on the Si surface. The energy resolution for this 
model is 3 eV FWHM, with a 10 eV threshold. Pulses from 
this model are shown in Fig 2] 



Parameter 


Value 


Units 


Description 


Csi 


43.3 


pJ/K 


Debye heat capacity 


Ctes 


31.1 


pJ/K 


TES electron heat capacity 


Cr cp 


29.3 


nW/K 


TES-Si thermal conductance 


G p b 


0.17 


nW/K 


Si-bath thermal conductance 


T b 


7.5 


mK 


Cold bath temperature 


T c 


15 


mK 


TES temperature 


R 


3 


mf2 


Quiescent TES resistance 


Io 


14.1 


fiA 


Quiescent TES current 


Po 


0.6 


pW 


Quiescent TES power 


R dT 


50 


- 


TES sensitivity 


To 


436.2 


ms 


Natural decay time Ctot/G p b 


T e a 


51.1 


ms 


Response time with TES speedup 


Tdccay 


29.2 


ms 


Decay time with readout circuit 


L 


30 


^H 


Readout inductance 



Multiplexing readout schemes for transition-edge sen- 
sors are now a mature technology being developed for 
many astronomical applications, for example 38|, and 
10,000 channel systems with time constants similar to 
this application are already in operation [391 ] . Schemes 
for even larger multiplexing gains arc in development [40J . 
Given the slow time constants of this application, a 
10,000 channel multiplexer design carries a fairly low risk 
and would allow 500 kg of Si to be instrumented. 

A concept for a 500 kg payload is shown in Fig. [5] 
The 10,000 Si bolometers are arranged in a column of 
dimensions 0.42 (dia.) x 2.0 (length) meters inside a di- 
lution refrigerator suspended from a vibration isolation 
mount. Passive or active shielding surrounds the refrig- 
erator. The exact shape and type of shielding will be 
determined at a later time. A cylindrical bore, perhaps 
10 cm or less in diameter, is removed from the shield and 
allows the 37 Ar source, mounted on a radio-pure transla- 
tion mechanism, to be moved to different positions along 
the side of the array. Periodic movement of the source 
throughout the measurement sequence allows each detec- 
tor to sample multiple baselines, enables cross-calibration 
among detectors, and aids in background subtraction. 
The minimum distance from the source to a bolometer is 
assumed to be ^10 cm. 
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FIG. 5: Conceptual schematic of the experimental setup 
for a bolometric measurement of coherent scattering from 
a high-intensity 37 Ar neutrino source. An array of 10,000 
Si bolometers is arranged in a column of dimensions 0.42 
(dia.) x 2.0 (length) meters (shown in green) inside a dilu- 
tion refrigerator suspended from a vibration isolation mount. 
Each Si bolometer has a mass of 50 g for a total active mass of 
500 kg. Appropriate passive or active shielding surrounds the 
refrigerator. A cylindrical bore in the shield allows the 37 Ar 
source, mounted on a translation mechanism, to be moved to 
different positions along the side of the array. Periodic move- 
ment of the source throughout the measurement sequence al- 
lows each detector to sample multiple baselines, enables cross- 
calibration among detectors, and aids in background subtrac- 
tion. The minimum distance from the source to a bolometer 
is assumed to be ~10 cm. 



Detector Backgrounds 



The detectors described in the previous section will be 
sensitive to several sources of background in the recoil 
energy range of 10-50 eV. Unfortunately, it is difficult to 
estimate accurately the rate of events from each of these 
sources, and the levels expected are currently unknown. 
Although we do not have a quantitative understanding 
of the backgrounds in this regime, it is important to note 
that backgrounds can be measured and subtracted us- 
ing data taken when the neutrino source is not in place. 
We qualitatively consider several sources of background 
which we expect to be present in the energy range 10- 
50 cV: 



• Radiogenic impurities: Most radioactive decay 
products have energies in the range of hundreds of 
keV to tens of MeV and will be clearly distinguish- 
able from the neutrino signal. Radiogenic impuri- 
ties may still contribute to the background in two 
primary ways. First, many common impurities pro- 
duce gamma rays that can interact with material by 
the photoelectric effect or Compton scattering and 
produce background events by the mechanisms de- 
scribed below. These gammas commonly arise from 
the U and Th chains and also from 40 K and 60 Co. 
Second, electrons from beta decay isotopes, such 
as tritium, may have arbitrarily small energies and 
therefore can produce electron recoils in the signal 
region. Nuclear recoils from decays at the detector 
surface, in which the electron is undetected, may 
also deposit small amounts of energy. 

• Compton scattering: Photons from radioactivity 
and atomic transitions in the detector material or 
housing may Compton scatter once at a shallow 
angle in the detector. Since there is no discrimi- 
nation between electronic and nuclear recoils, such 
shallow scattering would be indistinguishable from 
the neutrino signal. While the rate of these events 
is obviously dependent on the level of radioactive 
contamination, we expect kinematics to strongly 
suppress the rate of these events. 

• Photoelectrons: Photons produced in the detector 
or housing may produce photoelectrons in the de- 
tector material, which could be ejected. Recoils 
from such events could produce small energy depo- 
sitions in the energy region of interest. Low-energy 
secondaries from high-energy gammas produced in 
the detector or housing may impinge on other in- 
active material in the experiment and eject low- 
energy photoelectrons that could strike a detector. 

• Photons from atomic relaxation transitions from 
the surrounding surfaces: Photons from atomic 
transitions are of roughly the correct energy to pro- 
duce some background events near a 10 eV thresh- 
old in a Si detector. Copper, a good material for the 
detector housing for example, has 250 atomic lines 
with energies in the range 10-50 eV. The rate ex- 
pected from such events is very difficult to quantify, 
and would depend on the amount of low-energy ra- 
diation present in the cryostat to excite these tran- 
sitions. 

• Neutrons: Conventional methods used to reduce 
and model the neutron flux in dark matter ex- 
periments can be used. Running the experiment 
at large overburden, a muon veto can be used to 
veto cosmogenic neutrons with high efficiency The 
background rate of neutrons from muons that miss 



the veto can also be estimated with simulation. 
CDMS, for example, is able to achieve an unvetoed 
neutron rate of < 0.1 events/kg/year in the energy 
range lOkeV to lOOkeV [28|. Since the cross section 
for elastic scattering of neutrons on Si is fairly con- 
stant down to low energies, we do not expect the 
neutron background to be significant in our energy 
region of interest. 

• Neutrino- electron scattering: In addition to scat- 
tering coherently off nuclei, the neutrinos will also 
scatter off electrons in the detector material. The 
cross section for this process is lower than the cross 
section for coherent neutrino scattering, and the 
recoil spectrum extends up to hundreds of keV for 
37 Ar neutrinos. Since this background is also well- 
predicted by the standard model, we expect this to 
be a small contribution that may be reliably sub- 
tracted. 

• Dark matter: Recent experiments have reported 
signals that are consistent with a light WIMP of 
mass ~ 7 GeV and spin-independent cross section 
of ~ 10 -40 cm 2 for elastic scattering 4l|,|42j. If the 
dark matter interpretation of these signals is cor- 
rect, then scattering of WIMPs in the bolometers 
could be significant in the energy range of interest 
for coherent neutrino scattering. With standard 
assumptions regarding the WIMP halo velocity, we 
estimate a rate of 0.16 events/kg/day in 10-50 eV 
for such a - 7 GeV WIMP. Although this is a sub- 
stantial rate, it is much smaller than our expected 
event rate and would be independent of the dis- 
tance of each detector from the neutrino source. 
Given the uncertainty surrounding the measure- 
ments in [4l|, |42| and the conflicting result reported 
in [43[ , we do not consider this background further 
in our analysis. 

Excluding dark matter and the unknown backgrounds 
due to atomic transitions, the Compton scattering and 
photoelectron backgrounds are expected to dominate. 
Using the raw rate of events in CDMS, we can conser- 
vatively estimate the background rate due to these two 
sources. A good detector (250 g) in CDMS sees a raw 
rate of 0.001 Hz between 1-200 keV. If we conservatively 
assume that these events all lie in the range 1-10 keV 
and furthermore are attributable only to Compton scat- 
tering and photoelectric effect, then the rate of events 
is 38.4 events/kg/day /keV, assuming a flat spectrum. 
Assuming that the spectrum is flat down to 10 eV, we 
would see 1.54 events/kg/day in the energy range 10- 
50 eV. In reality, we do not expect the spectrum to be 
flat down to low energies. Low-energy gammas from k- 
and 1-shell electron captures are much more likely to be 
absorbed by the photoelectric effect, causing their full 



energy to be measured in the detector, and suppress- 
ing their background at lower energies. Only a small 
fraction of these produce ejected photoclectrons which 
could leave small energy depositions in the detector as 
described above. Compton scattering is very suppressed 
in the 1-10 keV range, and only dominates above 60 keV. 
The fraction of the background seen in CDMS detectors 
that lies above 60 keV is very small. For these reasons, 
the figure of 1.54 events/kg/day between 10-50 eV is a 
very conservative estimate of the Compton and photo- 
electron background. For the sake of simplicity, we as- 
sume a 1 events/kg/day background in our energy win- 
dow of interest. 



Detector Calibration 

Detector-to-detector variations constitute the most se- 
rious uncertainty in this measurement. Differences in the 
yield at different distance could quickly mask -or worse 
yield a false-positive on- the oscillation signal. Part of 
the problem can be solved by calibrating the efficiency of 
the detectors using a neutron calibration source. An at- 
tractive candidate would be to use a 3 H( 3 H,2n) 4 He com- 
pact di-neutron source, such as used in the petroleum 
industry [44|. The neutrons produced from this source 
should yield a relatively flat neutron energy spectrum, 
which is ideal for studying detector acceptance and to 
verify the 1/r 2 response of the array. 

As shown in FigureEJ our approach is to remove the de- 
pendency of the oscillation measurement on the response 
of a particular detector. The source is placed on a mov- 
able platform, and moved along the Si array through- 
out the measurement. Over the course of the experi- 
ment, each detector samples multiple baselines, and can 
be cross-calibrated with other detectors at each baseline 
to remove the individual detector response differences. 
Detector variations are essentially constrained by the in- 
situ measurement. 

It will also be important to calibrate the detector re- 
sponse to low energy photons and electrons. Recent ad- 
vances in solid state UV diodes make well-tuned eV pho- 
ton sources readily available. UV diodes ranging from 
255 nm to 350 ran (3.5 eV - 5.2 eV) with sub-eV resolu- 
tion are now commercially available. 



SENSITIVITY AND OUTLOOK 

Having discussed in detail both the source and the de- 
tection mechanism, we can now examine the signal in 
such an apparatus. For a monochromatic and isotropic 
source with activity R{t) encapsulated in some volume 
Vs, the signal rate as a function of time t is given by the 
expression: 



S(t)=J2R(t)-ME,)-f(E v ,T ).^.p t J^- J 



r 



P(E„,r st ) 
^r 2 st 
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(8) 



where Na is Avogadro's number, pt is the target density, 
dVr,i is the differential target volume of a single detector, 
and r s t is the source-target distance. The sum is taken 
over all discrete detectors available for the measurement. 
In the approximation of a point source, Equation [8] re- 
verts to the more familiar form: 



S(t) =Y,m ■ ME.) ■ f(E v ,T ) ■ ^ • M t P{Ev '' 



Airrf 



(9) 
where r^ now is the average-weighted distance from the 
source to the individual detectors and M t is the mass of 
each detector. In the limit where the measurement time 
is much greater than the source half-life, the total number 
of accumulated signal events is given by N ~ S(t ) ■ 
Ti/ln2, where r is the half-life of the neutrino source. 

2 ' 

Extending measurements well beyond the peak source 
activity has the added benefit of reducing the statistical 
uncertainty on the background. 

For a monoenergetic source, the oscillation signal is all 
encoded within the spatial distribution of events. A devi- 
ation from the expected r~ 2 dependence could constitute 
a possible oscillation signal. For the case where there is 
only one additional neutrino, the oscillation probability 
is given by the neutrino oscillation formula: 

P(E v ,r) = 1 - sin 2 (26 s) sin 2 (1.27 Am 2 s ^—) (10) 

E v 

where sin 2 (26$) is the amplitude to oscillate to the sterile 
state, and Am| represents the sterile mass splitting. In 
this case, E v is measured in units of MeV, r in meters, 
and Arrig in eV 2 . For simplicity, we look at the simple 
3+1 neutrino model, where the oscillation is to just one 
additional sterile neutrino. 

We use simulated data from a mock experiment to de- 
termine the potential sensitivity to sterile neutrinos. We 
consider a compact 5 MCi 37 Ar source to be used in con- 
junction with a 500 kg silicon array. We consider a total 
exposure of 300 days in order to extract both signal and 
background rates. Parameters relevant for the fit are 
listed in Table IIIII For comparison, we also list the pa- 
rameters for a germanium array with similar number of 
detectors and energy threshold. Due to the lower mass 
per detector needed to achieve the lower threshold and 
the lower recoil energies of neutrinos off the heavier ger- 
manium nucleus, a germanium array will achieve a signal 
rate that is about half of the silicon array. 

For such an experiment, we also consider a number of 
systematic errors: 



TABLE 111: List of relevant source and detector parameters 
used for sensitivity analysis. The signal rate is quoted for a 
single detector located 10 cm away from the center of a 5 MCi 
(185 PBq) 37 Ar source. 



Parameter 


Detector Type 


Detector Material 


Si Ge 


Atomic Number 


28 72.6 


ao(Eu) (KT 42 cm 2 ) 


0.44 3.82 


J- max 


50.3 eV 19.4 eV 


Threshold 


10 eV 


f(E v ,T ) (seeEq. @) 


64.2% 23.6% 


Detector cube size 


28 mm 15.5 mm 


Detector Mass 


50 g 20 g 


Number of Detectors 


10,000 


Total Mass 


500 kg 200 kg 


Yield at 10 cm (kg - May -1 MCi -1 ) 


15.28 19.0 


Signal Rate at 10 cm 


3.82 day" 1 1.90 day" 1 



• Source Strength: The SAGE collaboration used a 
variety of techniques in order to determine the fi- 
nal Ar activity, including gas volume, gas mass, 
calorimctry, direct counting and isotopic dilution. 
Any one of these methods in isolation achieved a 
±1% accuracy, while in conjunction the total un- 
certainty was less than ±0.5%. In this study, we as- 
sume a conservative ±1% uncertainty on the source 
strength. Since the source uncertainty applies to 
all detectors globally, it has minimal impact on the 
oscillometry measurement. 

• Cross-section: The cross-section uncertainty, much 
like the source strength uncertainty, is a global un- 
certainty and has little impact on our oscillome- 
try extraction. Its uncertainty would nominally be 
dominated by the uncertainty in the form factor, 
but at such exchange momenta the effect is ex- 
pected to be small. We therefore assume a ±1% 
global uncertainty due to the cross-section. 

• Vertex Resolution: The bolometric detector in this 
experimental design is composed of 10,000 silicon 
or germanium absorbers instrumented with a sin- 
gle thermometer. The dimensions of these absorber 
cubes are 28 and 15.5 mm per side for silicon and 
germanium, respectively. These dimensions arc 
smaller than the source itself (assumed to have a 
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radius of 4 cm) , thus the vertex resolution is dom- 
inated solely by the extension of the source. This 
effect is incorporated into our analysis. 

• Detector Variations: Detector variations are kept 
under control via the scries of in-situ and ex-situ 
calibration measurements discussed in the previ- 
ous section. Using the movable source depicted in 
Fig. El one should be able to calibrate the detector 
variations to about ±2%. The global uncertainty, 
which also depends on fiducial volume dependence, 
overall efficiency, etc., is estimated to be ±5%. 

• Detector Backgrounds: With detector-to-detector 
variations calibrated away, the main challenge for 
such a measurement remains the number of detec- 
tor backgrounds that accumulate during the mea- 
surement. As discussed above, the question of what 
the background will be between 10-50 eV is hard 
to estimate at this point, and more work needs to 
be done to enable a credible estimate. For this 
study, we assume a total background activity of 1 
event/kg/day in the signal region of interest. The 
signal-to-noisc ratio should scale roughly as the 
square root of the number of background events. 
The dependence of the accuracy of the measure- 
ment as a function of then signal-to-noise ratio is 
shown in Fig. HI Measurements taken with back- 
ground levels below 1 event/kg/day are essentially 
systematics dominated. 

• Source- Induced Backgrounds: Any backgrounds 
that stem directly from the 37 Ar source may po- 
tentially dilute the sensitivity of the measurement, 
since they, too, would exhibit a 1/r 2 behavior. The 
SAGE group has produced an extremely pure argon 
source, with less than 1% of the volume having 39 Ar 
contamination. Given the extremely long half-life 
of 39 Ar and the vastly different signature (/3-decay), 
we believe this is a negligible background source. 
Though the majority of the energy from the decay 
of 37 Ar is removed by neutrinos, a fraction of the 
energy is carried away from recoils and internal- 
brcmsstrahhmg photons. The SAGE source effec- 
tively reduced this contribution to less than 0.2%. 
Effective shielding should reduce this contribution 
even further. 

• Other Neutrino Interactions: One of the isotopes 
of germanium ( 71 Ge) has a low enough threshold 
to allow v e charged current scattering. However, 
with a threshold energy of 321 keV, the outgoing 
electron will have recoil energy far above the en- 
ergy region of interest. Hence, its contribution to 
the overall background can be considered negligi- 
ble. Charged current interactions on silicon all have 
thresholds above 1 MeV, hence they do not con- 
tribute to the background activity. As discussed 



previously, other charged-current interactions, such 
v e e~ elastic scattering, arc highly suppressed. As 
such, their contribution is also expected to be neg- 
ligible. 
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FIG. 6: Plot of the relative signal error versus signal-to- noise 
ratio -§— (S represents signal strength, B represents back- 
ground counts) for a 500-kg Si array exposed to a 5 MCi 37 Ar 
for 300 days. This array configuration and source intensity 
yields approximately S ~ 54, 000 total signal events. Arrow 
indicates signal-to-noise ratio corresponding to 1 background 
event/kg/day. 



TABLE IV: List of systematic uncertainties expected for a 
low-threshold germanium detector array. Uncertainties are 
listed for both shape+rate and shape only analysis. 



Source 


Systematic 




Global 


Shape Only 


Source Strength 


±1% 


- 


Cross-section 


±1% 


- 


Detector Variation 


±2% 


±2% 


Absolute Efficiency 


±5% 


- 


Source-Induced Background 


<1% 


<1% 


Vertex Resolution 


±2.8 cm 


±2.8 cm 


Source Extent 


±4 cm 


±4 cm 


Total Systematic 


±5.5% 


±2% 


Statistical (Whole Array) 


±1% 



A summary of the relevant systematic uncertainties 
are listed in Table Hvl A simple x 2 -fit is used to esti- 
mate the sensitivity of the proposed Si and Ge arrays. 
The data extracted from the entire array is first fit as 
a function of time in order to extract the overall source 
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strength and background (see, for example, Figure [7]). 
The background-subtracted signal is then fit to the oscil- 
lation formula of Eq. EH The analysis uses both shape 
and rate to determine the sensitivity to sterile neutri- 
nos. In the case of the shape+rate analysis, an additional 
penalty term is added to the likelihood from the overall 
flux measurement. 

Results for a 500 kg Si detector array are shown in 
Fig. 9(a) The distortion caused by a non-zero sterile 



mixing is statistically distinguishable in the measured 
distance profile (see Figure 8(a)). As can also be seen 



from the figure, the array is not necessarily fully opti- 
mized for a given oscillation length scale. Such optimiza- 
tion can proceed once the parameter space for sterile 
neutrinos is further constrained by ongoing and future 
neutrino experiments. Nevertheless, for the bulk of the 
region of Am| = 1-10 eV 2 and sin (20 s ) 2 > 0.08, typi- 
cally preferred from the reactor data, is ruled out at the 
90% C.L. If the best fit solution from the reactor anomaly 
is viable, then the measurement should be detectable at 
the 99% C.L. (sec Fig. |8(b)"| . It is possible to also conduct 
a shape-only analysis. Most of the sensitivity to sterile 
oscillations is retained for Ato| masses below 10 eV 2 . A 
sensitivity curve for a shape only analysis is shown in 
Figure EH 

For comparison, we also consider an equivalent Ge ar- 
ray with a total mass of 200 kg. These results are shown 
in Fig. |9(b)[ Finally, for completeness we also show the 
detector sensitivity for the Ge and Si arrays using an 
equivalent 51 Cr radioactive neutrino source (Figures 9(c) 



and 9(d) |. In general, the reduced source energy de- 



creases the available statistics, so a relatively stronger 
source needs to be considered in such a case. 
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FIG. 7: Distribution of events on a 500 kg Si array as a func- 
tion of time of source deployment. Source considered here is 
a 5 MCi 37 Ar electron-capture source. 



SUMMARY 

We have outlined the possibility of probing the exis- 
tence of sterile neutrinos using coherent scattering on 
a bolometric array. Such a method could provide the 
most direct test of oscillations to sterile neutrinos. With 
the advent of low threshold detectors and the use of in- 
tense neutrino sources, such an experiment appears fea- 
sible with our current technology. Such a program is also 
very complimentary to any existing dark matter search. 

Even in the absence of sterile neutrinos, the experi- 
ment as described in this letter can make other impor- 
tant measurements. Most prominently, such an exper- 
iment may constitute the first observation of coherent 
scattering. For a 500 kg detector, it should be able to 
make a ~ 5% measurement on the overall cross-section, 
pending on the absolute calibration of the efficiency. For 
an isoscalar target, such as silicon, this provides a direct 
measurement of the weak mixing angle at momentum 
transfer as low as 1 MeV. 
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FIG. 9: Likelihood contours for a 300-day run on a 500 kg Si array (left) and 200 kg Ge (right) array exposed to a 5 MCi 37 Ar 
(top) and 51 Cr (bottom) source, using both shape and rate information. Confidence levels in all plots are shown at 90% (blue), 
95% (green), and 99% (red). Statistical and systematic errors are included in the signal analysis. 
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FIG. 10: Likelihood contours for a 300-day run on a 500 kg Si 
array and a 5 MCi 37 Ar source, using only shape information. 
Contour curves use same scheme as in Figure [8] 



